Quantum Ice : a quantum Monte Carlo study 
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Ice states, in which frustrated interactions lead to a macroscopic ground-state degeneracy, occur 
in water ice, in problems of frustrated charge order on the pyrochlore lattice, and in the family 
of rare-earth magnets collectively known as spin ice. Of particular interest at the moment are 
"quantum spin ice" materials, where large quantum fluctuations may permit tunnelling between a 
macroscopic number of different classical ground states. Here we use zero-temperature quantum 
Monte Carlo simulations to show how such tunnelling can lift the degeneracy of a spin or charge ice, 
stabilising a unique "quantum ice" ground state — a quantum liquid with excitations described by 
the Maxwell action of 3+1-dimensional quantum electrodynamics. We further identify a competing 
ordered "squiggle" state, and show how both squiggle and quantum ice states might be distinguished 
in neutron scattering experiments on a spin ice material. 

PACS numbers: 74.20.Mn, 11.15.Ha, 71.10.Hf, 75.10.Jm 
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Ice is one of the strangest substances known to man. 
In the common forms of water ice, protons occupy the 
space between tetrahedrally coordinated oxygen ions, 
and each oxygen obeys the "ice rule" constraint of form- 
ing two long and two short bonds with neighbouring pro- 
tons P, [|| . It was quickly realised that these ice rules did 
not select a single, unique proton configuration but 
rather a vast manifold of classical ground states, with an 
entropy per water molecule of s ~ ke log(3/2) @. This 
prediction proved to be in good agreement with measure- 
ments of entropy at low temperatures Q, but stands in 
clear violation of the third law of thermodynamics — at 
zero temperature we expect water ice to be described by 
a single, unique, ground state wave function. 

The same "ice" rules, and the same extensive ground 
state degeneracy, arise in (i) problems of frustrated 
charge 0, @ and orbital @ order; (ii) proton-bonded 
ferroelectrics [3]; (in) statistical descriptions of polymer 
melts Q; and (iv) a family of rare-earth magnets col- 
lectively known as "spin-ice" 0-12|. In each case, the 
ice rules have non-trivial consequences for the proper- 
ties of the system, notably an algebraic decay of correla- 
tions |71 |13h15|| and excitations with "fractional" charac- 
ter fllM la. I la [13] ■ These exotic features of the ice state 
have been extremely well characterized in spin ice, where 
the algebraic decay of correlation functions is visible as 
"pinch points" in the magnetic structure factor 18| , and 
the fractional excitations have the character of magnetic 
monopoles (ljil - t^ . 



All of these systems beg the obvious question — how 
is the degeneracy of the ice manifold lifted at zero tem- 
perature ? The simplest way for an ice to recover a 
unique ground state at zero temperature is for it to or- 
der. This is exactly what happens in KOH-doped water 
ice, where the protons order below 70K [23j]. However in 
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FIG. 1: (Color online) Quantum mechanics enters ice physics 
through the tunnelling of the system from one ice configuration 
to another. The leading tunnelling matrix element for spin ice is 
the reversal of a set of Ising spins with closed circulation on 6-link 
hexagonal plaqucttc. 



many spin-ice materials, no order is observed [121 ]. This 
raises the intriguing possibility that there might exist a 
zero temperature "quantum ice" state, in which a single 
quantum mechanical ground state is formed through the 
coherent superposition of an exponentially large num- 
ber of classical ice configurations. Such a state could 
have a vanishing entropy at zero temperature, and so 
satisfy the third law of thermodynamics, without sacri- 
ficing the algebraic correlations and fractional excitations 
(magnetic monopoles) associated with the degeneracy of 
the ice states. 

In this Letter we use zero-temperature quantum Monte 
Carlo simulations to establish the ground state of the 
minimal microscopic model of a charge or spin ice with 
tunnelling between different ice configurations. We find 
that the ground state is a quantum liquid, with an emer- 
gent £7(1) gauge symmetry, and excitations described by 
the Maxwell action of 3+1-dimensional quantum electro- 
dynamics. This state is the exact, quantum, analogue of 
the spin-liquid phase realised in "classical" spin ices such 
as Dy2Ti20"7, and exhibits the same magnetic monopole 
excitations. We also explore how quantum effects in this 
novel liquid modify the "pinch-point" singularities seen 
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in neutron scattering experiments on spin ices. 
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FIG. 2: (Color online) Ground state phase diagram of the 
quantum ice model [Eq. l(2j] as a function of the ratio fi of kinetic 
to potential energy. The 3D "quantum ice" point fi = exists 
deep within an extended quantum liquid phase with deconfincd 
fractional excitations and algebraic decay of correlation functions. 



The best systems in which to look for a quantum ice 
are those which are able to tunnel from one ice configu- 
ration to another. In water ice, in the absence of mobile 



ionic defects [24[ , this tunnelling occurs through the col- 
lective hopping of protons on a 6-link loop. In spin ice 
it is the c ycli c exchange of Ising spins on a hexagonal 
plaquette [251] . illustrated in Fig. LTJ In both cases the 
ice rules can be written as a compact lattice U(l)-gauge 
theory in which the displacement of protons — or orien- 
tation of magnetic moments — are associated with a fic- 
titious ma gne tic field B = V x A, in the Coulomb gauge 
V • A = [13j, [l5| . Tunnelling between different ice con- 
figurations introduces dynamics in the gauge field A, and 
the minimal description of a quantum ice is the Maxwell 
action of conventional quantum electromagnetism 



J d 3 xdt [E 2 - c 2 B 2 ] 



(1) 



where c is the effective speed of light and, in the absence 
of electric charges, E = —dA/dt. It follows directly from 
Eq. (JTJ that correlation functions have dipolar character, 
and local defects in an ice configuration act like decon- 
fined magnetic monopoles 26H28j. 

Field-theoretical arguments alone cannot resolve 
whether the quantum t/(l) liquid described by Eq. (JTJ 
is realised in an "ice" material, or realistic microscopic 
model. Encouragingly, evidence supporting the existence 
of such a phase has been found in quantum Monte Carlo 
simulations of strongly-interacting hard-core Bosons on a 
pyrochlorc lattice [29( . However, these simulations are re- 
stricted to temperatures comparable with the degeneracy 
temperature of the ice manifold, and so are mute as to 
the zero-temperature ground state. In this article, we use 
zero-temperature Green's function Monte Carlo (GFMC) 
simulation techniques [lo[ to provide concrete evidence 
for the existence of a quantum U(l) liquid ground state 
in a microscopic lattice model of a quantum ice. 

The model we consider was first introduced by Her- 
mele et al. 27] as an effective Hamiltonian for an easy- 




FIG. 3: (Color online) The ice configuration posscsing the most 
flippable plaquettcs is the "squiggle" state, shown here within a 
40-site tetragonal cell. Arrows show the displacement of protons 
within Ic water ice or, equivalently, the orientation of spins in spin 
ice. The squiggle state posses a net "magnetic" flux, <j> = <^ squ iggi e , 
orientated along a [100] axis of the crystal. In a spin ice material 
this would correspond to a net magnetisation of 1/5 the value at 
saturation, directed along a [100] axis. 
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acting on all possible (spin) ice configurations.. The first 
term in Eq. ^ describes tunnelling from one ice con- 
figuration to another, where | O) should be understood 
as a closed circulation of B on a "flippable" hexagonal 
plaquette [cf. Fig.Q]. The sum X/piaq runs over all such 
plaquettcs in the lattice. The additional potential en- 
ergy term /x counts the number of flippable plaquettes in 
a given ice configuration, and renders the model exactly 
soluble for /x 
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All energies are measured in units 
of the tunnelling matrix element between ice configura- 
tions. For /i = 0, Eq. ([2]) is the minimal microscopic 
model for a 3D quantum ice. 

We have previously used GFMC simulation to estab- 
lish the existence of a quantum U(l) liquid in the quan- 
tum dimer model on a diamond lattice [32|, [33|. The 
quantum ice model Eq. [2] differs from this only in that 
the Hamiltonian acts on fully-packed loop, rather than 
dimer coverings of the lattice [cf. Ref. |34| . We can there- 

We 
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axis antiferromagnet on a pyrochlore lattice. It is defined 



fore solve it using the methods set out in Ref. 
consider clusters with periodic boundary conditions, and 
make extensive use of the fact that the "magnetic" flux 
<f> = J dS ■ B through any periodic boundary is a con- 
served quantity. This makes it possible to define a series 
of flux quantum numbers <f> = {4> x , <p y , <p z ) for each clus- 
ter (33( 1 . Our findings are summarised in Fig. [2] 

For fx — > —00, the ground state of Eq. © is the ice 



3 




0.6 -0.4 -0.2 



FIG. 4: (Color online) Evidence for a phase transition out of 
the squiggle state, (a) Melting of squiggle order within the squig- 
gle flux sector <j> = squ iggie- The relative number of flippable pla- 
quettes v = Nf(fi)/Nf(—oo) is calculated using quantum Monte 
Carlo simulations for clusters of 160 to 1080 sites. The solid lines 
show different Pade approximants to the series expansion about 
a perfectly-ordered squiggle state, (b) Ground state level cross- 
ing between the flux sector associated with the squiggle state, and 
the zero-flux sector associated with the quantum U(l) liquid. The 
ground state energy per site is calculated for clusters of 432, 1024 
and 2000 sites (squigglc-flux sector) and 160, 540 and 1280 sites 
(zero-flux sector), using quantum Monte Carlo simulation. A clear 
crossing is observed for fi = —0.50 ± 0.03. 
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FIG. 5: (Color online) Evidence for the existence of a quan- 
tum U(l) liquid, (a) Flux dependence of finite-size energy gaps 
E^ — Eq for parameters bordering the RK point, fj, = 1. Results 
are plotted for quantum Monte Carlo simulation of an 1280-site 
cluster (solid points) and perturbation theory about the RK point 
(open black circles) . The dashed line indicates the scaling expected 
for a quantum U(l) liquid, following Eq. JT}- (b) Flux dependence 
of finite-size energy gaps at the quantum-ice point, \i — 0. Results 
are obtained using quantum Monte Carlo simulation for 320-site, 
640-sitc and 1280-site clusters, and in exact diagonalization for an 
80-site cluster. The dashed line indicates the scaling expected for 
a quantum (7(1) liquid. 



configuration which maximizes the number of flippable 
plaquettes. This is the 60-fold degenerate "squiggle" con- 
figuration shown in Fig. It is ordered, and therefore 
exhibits Bragg peaks in diffraction experiments j35j . A 
good measure of squiggle order is the relative density 
of flippable plaquettes v = Nf(fi)/Nf(—oo). In Fig. |U[a), 
we compare GFMC calculations of v for a series of finite- 
size clusters, carried out in the squiggle flux sector, with 
Pade approximants to a series expansion in 1/^ about 
perfect squiggle order. The agreement between the two 
calculations is essentially perfect for fi < —0.75, and the 
marked suppression in the number of flippable plaque- 
ttes for /j > —0.3 is strongly suggestive of the melting 
of squiggle order. However, perturbation theory about 
the exactly soluble "RK" point fi = 1 dictates that the 
ground state of Eq . (|p should be in the zero-flux sec- 
tor for [i — » 1 27J, |33j. And in fact the gentle demise 
of squiggle order is pre-empted by a ground state level 
crossing between the squiggle and zero-flux sectors at 
fi = -0.50 ± 0.03, shown in Fig. gjb). 

These results are consistent with a first order phase 
transition out of the squiggle state at fi = —0.5, but do 
not yet confirm the existence of the quantum U(l) liq- 
uid we are seeking. Fortunately Eq. ([!} also makes spe- 
cific predictions for how the quantum U(l) liquid phase 
evolves out of the RK point at fi = 1 [26|, r27l . |33| . Specif- 
ically, the finite-size energy gaps E$ — Eq should grow as 
E<j> — Eq ~ c 2 r/> 2 / L where E^ is the energy of the ground 
state with flux $ and c 2 oc + . . .fUHH. In Fig.^a), 



we present simulation results for E$ — Eq for fi < 1. 
We find good agreement between GFMC simulations and 
perturbation theory about the RK point 27, 33|, and a 
near-perfect collapse of both data sets according to the 
prediction of Eq. This confirms the existence of a 
quantum t/(l) liquid bordering the RK point, as pro- 
posed in [27| . 

However so far as real materials are concerned, the 
most interesting point in the parameter space is the 
"quantum ice" point fi = 0. Does this also conform to 
the behaviour expected of a quantum [/(l)-liquid ? In 
Fig. [S^b) we present simulation results for the finite-size 
energy gaps of Eq. @ at fi = 0. We extract the leading 
dependence on system size by plotting (E^ — E Q )/N s \ tes 
as a function of (0/0 sq uiggie) 2 - Once again, the collapse of 
the data is excellent, confirming the existence of a quan- 
tum J7(l)-liquid. 

The ground state phase diagram of the quantum ice 
model Eq. ^ is summarised in Fig. [21 For fi < —0.5, 
the ground state is the complex "squiggle" order shown 
in Fig. [31 For fx > —0.3 we find unambiguous evidence 
for a quantum liquid phase which is well-described by 
the U(l) (lattice) gauge theory of quantum electromag- 
netism, and so will exhibit both algebraic decay of cor 
relations and deconfined fractional excitations 
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This quantum U(l) liquid phase terminates in an RK 
point at fL = 1, and therefore the "quantum ice" point 
fi = lies deep within it. For —0.5 < fi < —0.3 sim- 
ulation results depend in detail on the geometry of the 
cluster chosen. However we tentatively conclude that a 
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single, first order phase transition takes place between the 
squiggle and J7(l)-liquid phases for u = —0.5. This phase 
diagram should be contrasted with those for bosonic [36| 
and fermionic j37j quantum ice models on the 2D square 
lattice, where all phases are ordered and confining for 
fi < 1. In short — 2D quantum ice models are ordered 
and confining, but the 3D quantum ice model solved here 
is not. 

So far as electronic charge ices are concerned, this work 
should be regarded as a "warm-up" exercise, since Eq. ([2]) 
does not allow for the spin or Fermi statistics of the elec- 
trons 33,37,33- Similarly, the application of these ideas 
to hexagonal water ice depends crucially on the general- 
isation to a different lattice, and the role of more general 
ionic defects [24[. However the model we have solved 
may give a good account of "quantum spin ice" mate- 
rials such as^bzTizOr |2§G3, Pr 2 Sn 2 7 0,0 and 
Yb 2 Ti 2 7 " 



l48j . where quantum fluctuations of 
magnetic moments provide a route to tunnelling between 
spin ice states. And in this context it is interesting to ask 
how a quantum ice might be distinguished in experiment 
on a spin ice material ? 
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FIG. 6: (Color online) (a) Static structure factor 5(q) in the 
spin-flip channel for a classical spin ice, as measured by Fennel et 
al. [HI , calculated here from the microscopic model Eq. J2j, with 
fj, = 1. Results are plotted in the (h,h,l) plane, and show the 
"pinch-point" structure associated with 1/r 3 dipolar correlations 
in 3 dimensions, (b) Equivalent static structure factor S(q) for a 
quantum spin ice described by Eq. JJJ with fi = 0. In this case, 
correlations between spins show the 1/r 4 behaviour characteristic 
of dipoles in 3 + 1 dimensions, and the pinch-points are eliminated. 
All simulations were performed for a cubic cluster with 2000 lattice 
sites. 



The signal feature of a classical ice state is the pres- 
ence of "pinch point" singularities in the static structure 
factor S(q) [IS]. These reflect the fact that the spins or 
charges which make up the ice state have correlations 
of 3-dimensional dipolar form 0, 0, EH. In contrast, 
in a quantum ice, static correlations take on the form 
of dipoles in 3 + 1 dimensions ana and as a result, 
the pinch-point singularities in S(q) arc eliminated. To 
illustrate this, in Fig. |6] we present GFMC simulation 
results for S(q) for a quantum spin ice, calculated di- 
rectly from Eq. $Z§. Results for the RK point fj, = 1, 
where the correlations are classical, clearly show pinch 
points at reciprocal lattice vectors. However at the quan- 
tum ice point, [i = 0, there is a marked suppression of 
spectral weight around the same reciprocal lattice vec- 
tors. As a consequence, the angle-integrated structure 
behaves as S(|q| —> 0) oc |q|, while in a classical spin ice 
S(\q\ — >• 0) —> const. This linear-|q| behaviour at small 
|q| is consistent with published results for Pr 2 Sn 2 C>7 [Z3I ]. 
It might also be interesting to reexamine elastic neutron 
scattering data on other pyrochlore antiferromagnets in 
the light of these results 1491 . These issues will be ex- 
plored further elsewhere [501 ] - 
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